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(1) The resonator-decay and the reverberation methods for measurement of the absorption of 
sound in liquids are in wide use, but their theoretical foundation has not hitherto been investigated. 
The basic assumptions are in fact invalid, but under most practical conditions the errors are tolerable, 
at least for three common geometries. (2) The very large excess losses always observed at the lower- 
order modes are identified with anelastic effects in the envelope, for the most part, and with viscous 
boundary-layer effects to a lesser, but very substantial extent. For modes of sufficiently high order 
these become negligible. 
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1. Introduction 

There is little doubt that many of the best measure- 
ments of the absorption of sound in liquids, especially 
at the lower frequencies, have been made by "decay" 
or "reverberation" methods. In the first case, a 
resonator formed by the sample liquid and a suitable 
container is excited in a single mode; the attenuation 
coefficient, a, is calculated from the decay rate after 
the excitation is cut off. The second case is much the 
same, except that the excitation is over a band of 
frequencies. 

It is curious that it has always been tacitly assumed 
that a as measured on a resonator is the same as that 
for a plane running wave in free space. That this cannot 
be true in general is evident from an obvious counter- 
example—in the case of uniform expansion (wave- 
length large re dimensions) the losses are governed 
by the congressional viscosity, \+2//,/3, whereas 
in free space the losses are governed by A+ 2/jl. 

In section 2 we analyze the decay for a single mode 
in general and apply the results to three common 
cases: radially symmetric vibrations of a sphere, 
first used by Leonard [1946] 1 and in wide use ever 
since; radially symmetric vibrations of a cylinder, 
apparently first suggested by Meyer and Tamm 
[Mudlers, 1948] and the general vibrations of a rec- 
tangular parallelepiped (box). It is found that for the 
geometries analyzed the assumption in question is 
in general false, although for suitable configurations 



1 Figures in brackets indicate the literature refe 



at the end of this paper. 



the errors in measured values of a may well be small. 
But even if so, it seems that a method in such wide 
use, and upon which so much reliance is placed, should 
rest on a firmer foundation than it has hitherto. 

Once precautions are taken to minimize radiation 
of sound from the resonator into the air and conduction 
of sound in the supports, the most important sources 
of error are probably dissipation in the envelope and 
dissipation in the boundary layer. Both Leonard [1950] 
and Wilson [1951] tried to calculate envelope losses 
in a sphere without much success. New calculations 
are presented in section 3. Boundary layer losses occur 
whenever the motion of the liquid near the wall has a 
component parallel to it. This was first pointed out to 
Leonard by R. J. Christensen, and it is on this account 
that radial modes of a sphere, for which losses are 
smaller than for any other configuration, were em- 
ployed by Leonard and by so many others since. 
Elementary calculations of boundary layer losses, 
for which no great originality is claimed, are presented 
in section 4 for comparison purposes. 

2. Theory 

The rationale of the resonator decay method can 
be summarized as follows. In a plane running wave, 
the energy density is 

E = E e- 2a * 

in which a is the amplitude attenuation coefficient. 
The variation with time, for an observer moving with 
the wave at speed c, is 

E = E e- 2act . 
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Logarithmic differentiation gives for the instantaneous 
fraction of E lost per unit time 



1 dE 



(1) 



This notion is now applied to a resonator. The total 
energy density is T + W, i.e., the sum of the kinetic 
and potential energy densities. The viscous energy 
loss per unit volume (heat conduction is here neglected 
for simplicity; for most liquids it is minor) is twice 
the dissipation function, 2F. The ratio — 2F/(T-\-W) 
corresponds to the left side of eq (1), hence 



a = 



J_ 2jFdV 1 jFdV 

2c ${T+W)dV~ 2c fWdV 



(2) 



the integrations being taken over the whole volume 
of liquid. The second expression on the right in eq 
(2) is not exact; its use will be justified shortly. The 
assumption we wish to test is that a in eq (2) is the 
same as that in eq (1). 

In the foregoing, and in the sequel as well, all 
first-order quantities, such as pressure, are supposed 
to vary as cos (ot. Second-order quantities, such as 
7\ W, and F, will have terms which vary as cos 2 o>£, 
sin 2 (L>£, or sin cot cos cot. We are interested in time 
averages of 7\ W, and F only — integration over t 
will yield coefficients of i or zero, as the cause may be. 
We bear this in mind, suppress the factors which 
indicate time variation, and take all first-order quan- 
tities as root-mean-square. 

The problem reduces to finding expressions for 
T, W, and F. These are gotten from the pressure 
and velocity fields, supposed known. Although these 
can be determined in principle (for simple shapes) 
the solutions are too cumbersome to be useful. In- 
stead we make the approximation that the field is 
the resultant of two others, the first appropriate to 
an inviscid fluid (i.e., one having a potential) and the 
second, appreciable only in a relatively thin boundary 
layer, that of a locally plane shear wave launched 
at the boundary along the inward normal, the magni- 
tudes of the two being adjusted so that the net tangen- 
tial velocity on the liquid-solid interface is zero 
[Morse and Ingard, 1968]. The functions T, W, and 
F are calculated from the potential-solution only. 

The assumed velocity potential, i//, satisfies the 
wave equation 



V 2 <// = -A: 2 i// = - (co 2 /c 2 )i//. 
The pressure is 

dip . , 



(3) 



(4) 



The comma indicates differentiation. We need also 
the strain rate 

Sij = i {iM,j + U4j)=—\lf,ij (6) 

and the dilatation rate 

A = Su = -ty,a = - V 2 i// = £ 2 i//. (7) 



The summation convention is used throughout. 

We are now in a position to calculate T 9 W, and 
F. First we note the constitutive equation 



Tij= (-p + kk)8 ij + 2iJLSi j 



(8) 



in which Ty is the stress, X and jjl the second and first 
viscosity coefficients respectively, and 8,j the unit 
tensor. The strain energy density is given by 

2W = TijSij = (- p + \A)SySy + 2/xS S/> (9) 

In eq (9) Sy is in quadrature with Sij and with A; the 
only surviving term is that in pbijSjj which equals 
pA, and as A = A/ico, we have from eqs (3) and (4) 



2W = pk 2 x\s 2 . 



(10) 



The kinetic energy density is pUiUi/2 = p^i^/2, from 
eq (5). A special case of Green's theorem, 

f^jfijijOr+f ^^dv=-[ ijjj^ds, (ii) 

leads to the condition of resonance. It is easily seen 
that it is equivalent to 

b T - W)dV =- P 2^fn dS - (12) 

At resonance, the total kinetic and potential energies 
are equal so that the right-hand side of eq (12) vanishes. 
A sufficient condition is that everywhere on the outer 
surface (of the container, which we treat as a fluid 
without too much error because it is relatively thin) 
either the pressure (proportional to i//) or the normal 
velocity (proportional to difj/dn) is zero. In practice, 
the resonator is operated at frequencies much above 
resonance (breathing mode) of the empty envelope. 
Hence the energy in the envelope is mostly kinetic 
and therefore most of the potential energy is in the 
liquid. A good approximation to the total energy is 
therefore 2/ v WdV, the integration being taken over 
the volume of the liquid. This explains the substitu- 
tion of 2JF for T+ W in eq (2). 

The dissipation function is given by 

2F=T ij S i j= (-p + \A)§ j S + 2/xS j Su- (13) 



in which p is the undisturbed density, and the velocity 
is 

Ui = — \jf,i (i = 1,2,3). (5) Here p is in quadrature withS^ and eq (13) reduces to 
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2F=(\ + 2M)A 2 + 4^(Sf 2 +S^+5|,-S„S 22 
- SiiSm - S-,;S r i ) = ( X + 2/ul ) AV + W-2 



(14) 



in which L> is the quadratic strain-rate invariant. 
Equation (14) together with eqs (2) and (10) gives 
finally, 

JFdV = (A + 2 M )o; 2 2^J/ 2 ^ 

2c/jjw 3>c 3 pck'frdv (15) 

We recognize the first term on the right-hand side of 
eq (15) as the plane-wave value of a; therefore the con- 
dition that the resonator-decay experiment yield the 
plane-wave value of a is that 

JhdV=0. (16) 

For reference we repeat 

/2=if 2 + S 2 2 3 + S3 2 1 -SnS22-S22S33--S33Sn. (17) 

We note two special cases. In the case of a plane 
wave, running or standing along ;n, only Sn in eq (17) 
is not zero, so I> = and eq (15) gives the known result. 
In the case of uniform expansion, the shear rates in 
eq (17) are zero and Sn =S 2 2 = S;w — S (say), uniform in 
space. Then A = 3S and I> = -35 2 = -*A 2 . With this 
substitution, eq (15) becomes 

(\ + 2/*/3)o> 2 
2pc* 

which is the correct result, \ + 2/x/3 being the compres- 
sional viscosity. 

2.1. Sphere, Radial Mode 

Here we use spherical coordinates, r, 0, (/>, and sup- 
pose that the potential, i//, depends on r only. Then the 
velocities are 



u r = - 



dtjj 

dr 



, uh — 11$ = 0. 



The nonzero strain rates are 



Srr = - 



d 2 l\) 



dr 1 
so that, from eq (17) 



1 difj 

,000— Ojxji — — - , 

r dr 



U r 



7 2 d±(d± + 2r d 



r 2 dr \dr 



dr 2 )' 



and with dV=^Trf z dr, 



fhdV=-*7T 



<f)' 



= — ^irau 2 , , 



(18) 



in which a is the radius of the sphere and u a = (u r ) r - 



Equations (18) and (15) together show that for this 
case the resonator-decay result underestimates a ex- 
cept only for u a = 0, corresponding to a rigid container. 

According to eq (14) for 2F the excess in a, Aa, can 
be estimated from 



Aa 
a 



ifiSM v 



(X + 2ix.)k*S v ^dV 



(19) 



For the radial mode of a sphere the velocity potential 
for unit amplitude is 



* = 



sin kr 



(20) 



Equations (18), (19) and (20) yield 



Aa 
a 



6 OLc\ 



ka — 



sin ka 
ka 



r 



ka a ka — sin ka cos ka 



(21) 



In eq (21) we have written a c \/a for 4/x/3(\ + 2/jl) 
where a c \ corresponds to the "classical" relationship 
3A-h2/x = 0. To use eq (21) for a particular case we 
need the values of k corresponding to the various 
resonances. Although these are readily calculated, 
we content ourselves with an examination of the worst 
case, i.e., the pressure-release case, p (l = Q. Then in 
eq (21) sin ka = 0, ka = mr, cos ka = ±\. We have 



Aa = _ ^ Ctrl 



I 



a n c TT- 



(22) 



In the case of water, for which a c \/a is about i, the 
"error" is large for the very low-order modes; about 
20 percent for n=l, five percent for n = 2, etc., but 
is less than one percent for n ^ 5. For most other 
liquids Aa/a is even smaller. As we shall see, there 
are other, more compelling, reasons to operate at the 
high-order modes, so that for practical purposes the 
spherical resonator, other things aside, has a reason- 
ably firm theoretical basis. 

As a check on eq (21) we note that for A: = (uniform 
expansion) eq (21) yields the known result. 2 

2.2. Cylinder, Radial-Axial Mode 

This case is not so simple, and the explicit expres- 
sion for if; is required. If the mid-point of the axis is 
the origin of cylindrical coordinates, r, 6, z, the 
potential is, for amplitude unity, 



*p = Jo(k f r) cos (kzZ + e), 



(23) 



with the radial and axial wave numbers k r and k z 
subject to 



k 2 + k 2 = k 2 = co 2 lc 2 . 



(24) 



2 The indeterminate form in eq (21) is easily evaluated by expansion of the sine and 
osine functions about ka = Q. 
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The velocities are 

U r = — — = krJ\C, Ue = 0, Ug = — — = kjos. 
dr oz 

Here and in the sequel we use abbreviations like J\ 
for Ji(k r r) and c for cos (k z z + e). The strain rates are 



Srr=^ = t*(j0-£)C, 



C Ur 1 2 J 1 

oee= — = £; — c, 

T KrT 



S zz = — 1 = k 2 JoS, 
oz 

• 1 (bU e dU r \ 

S r d = Sdz == 0. 

These, together with dV = 2irrdrdz and the limits to 
a for r and — L to L for z, give after some reduction 

--T- I WF=i<+fl [sin (fel + c) cos {k z L + e) 

TTkr J 

— sin (— k z L + e) cos (— /c^L-he)] 



where 



and 



A = 2LJ l (k r a) \k 2 z aj (k r a) 

B = k z a{k r aUS(k r a)+J 2 (k r a)] 
-Mk r a)Mkra)} + -^J*(kra) 



(25) 



drical surface [/i (A: r a) = 0] for which eq (26) becomes 

(27) 



a 



K,m 



a [ 2 ( niraY 



in which Xo,m is the tfith zero ofjo(^). Thus high mode 
numbers m or n (or both) insure a small "error." The 
conclusion is the same as for the sphere — the method 
is valid in practical cases. 

2.3. Rectangular Tank 

In rectangular coordinates, %\, the potential, for 
amplitude unity, is 



\jj = CiC 2 c 3 . 



(28) 



Again, we use abbreviations like C\ = cos (k\X\ + €1). 
In eq (28) the k\ are subject to 



We have 



Ui = 



and 



k\ + k\ + kl = k 2 = (w/c) 2 . 



difj 



= k 1 s 1 c 2 c 3 , u 2 = — - — = &2C1S2C3, 

OX\ 6X2 



u 3 = — ~ — =A:3CiC2S3, 
dx 3 



S n = — = A: 2 CiC 2 c 3 , etc. 
• 1 (du 2 . duA 

2 la — a/ = kxk2SxS2C ^ etc * 

The invariant / 2 of eq (17) will have three terms. The 
first is 

— - k\k\c\ [cos 2{k\x x + ei) + cos 2(fe*2 + e 2 )]; 



For the integral in eq (25) to be zero it is sufficient 
that Ji(k r a)=0, that either sin (k z L -he) or cos 
(k z L + e) be zero, and that either sin (— k z L + e) or 
cos(— k z L + e) be zero. This follows from the inde- 
pendence of k r and k z . Restated, a sufficient condi- 
tion that the experiment measure the free-space value 
of a is that the cylindrical wall be rigid and that each . ., , . (hJ . 

end be either rigid or pressure-release. (Another con- sin W Li + 6i) cos ^ + €i) 
dition is /c r =0, corresponding to a plane wave.) 

We can calculate the error, Aa/a, in the same way 
as for the sphere. In the special case of pressure- 
release ends, we get 



the other two are got by cyclic permutation of the 
subscripts. Integration over the range —Li to Li for 
Xi yields an expression for JhdV containing terms of 
the form 



— sin (— kiLi + €i) cos (— kiLi + €j) 



Ac* _ OLc\ k r 

a a k 4 a 2 



Jiikra) 



k z aj (k r a) + -£ Ji(k r a) 



together with cross-products of these taken two at a 
time. The coefficients of these are various products 
of Lj and powers of ku An account of the independence 
of these coefficients a sufficient condition that 



JZ(k r a)+J?(k r a) 



(26) 



fvh = 



is that each face be either rigid (sine term=0) or 
Again, the worst case is for a pressure-release cylin- pressure-release (cosine term = 0). There are, of 
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course, other special conditions. For instance if any 
two kt are zero the integrand I? vanishes — we are back 
to a plane standing wave. 

The actual "errors" for any particular case could be 
calculated with some labor, but it seems hardly worth- 
while in consideration of what we have already learned 
in the other cases. 



3. Losses in the Envelope 

3.1. Sphere, Radial Modes 

The loss in a spherical shell oscillating in the radial 
(breathing) mode is not difficult to calculate if a simple 
expression is obtained for the real component of the 
impedance of the shell. It is the fraction \/Qe of the 
stiffness reactance (not the mass reactance; the losses 
are associated with strain rates, not with inertia). The 
stiffness reactance is M^coo/ft), Me being the mass of the 
envelope and con its resonance angular frequency when 
empty. Hence the dissipation rate is 



1 — ka cot ka = 



piA 



D E = 



M E a>l 

Qe(*> 



(29) 



in which u (l is the radial velocity at r= a, the interface. 
The approximation is here made that the radial velocity 
is constant throughout the thickness of the envelope. 
This expression is of use only if the variation of Qe with 
co is known. For many solid materials, including glass, 
the Q at frequencies of interest here is independent of 
frequency (elastic hysteresis) [Mason, 1950J and we 
shall take it as such. We seek the ratio of De, the dissi- 
pation rate in the envelope, to De; that in the liquid. 
For the latter we use, without significant error only the 
part (k + 2iJL)k 4 f v ip 2 dV [eq (14)]. This and the value of 
u a in eq (29) are gotten from the potential in eq (20). The 
result is 



De _ 2p E CLt sin 2 ka /coo\ 2 (1 — ka cot ka) 1 
Dl (X-f2/x)(^;(A:a) 2 \ co/ ka — sin ka cos ka ' 



(30) 



in which ps and t are the density and thickness of the 
envelope and ci is the speed of sound in the liquid. 
Equation (30) is meaningful in our context only for those 
values of co (and k=o)/c) which correspond to actual 
resonances. 

One expression of the condition of resonance is that 
on the boundary between the liquid and the envelope 
the point admittance of the liquid equals that of the 
envelope. The former is 

Uu _ — c)ifj/da _\ — ka cot ka 



pa i(x>pi)\f (l i(x)piA 

Pl being the density of the liquid, and the latter is 



y E = llmp E t 



1 



coo\- 

CO 



(31) 



PEt 



(32) 



This together with eq (30) gives, after much reduction 



1 



(Oq 



De 



D L ceQe<*> ( ot-l co- ) \ CO 



R 



/? 2 -fl[l-(Wco) 2 ] + (M 2 [l-(Wco) 2 ] 2 ' 



(33) 



in which we have written R for pEa/p E t, and have elim- 
inated X + 2^t in favor of ct/co 2 from eq (15). We shall see 
that in most practical cases co <^ co so that the brack- 
eted expressions in eq (33) are nearly unity. Equation 
(33) is still not free of co (or k=co/c), but the right- 
hand side is seen to be a monotonic function of co, so 
that the equation is usable as it stands; it is not nec- 
essary to solve the frequency equation, eq (31) 
explicitly. 

It seems clear that, everything else being equal, 
there is no optimum thickness, t, of the envelope wall, 
because the wall loss would vanish at both extremes — 
no wall and rigid wall. This notation is verified by eq 
(33) according to which D e /De has a maximum for 
R = ka[l — (coo/co) 2 ]. (R varies as 1/0- The variation 
of the error, D E /Di with frequency is very rapid. For 
a given set-up, and for liquids for which cx/co 2 is 
independent of frequency, eq (33) predicts that for 
low-order modes, such that R > ka, De/De varies as 
l/co\ and for high-order modes, such that R < ka, it 
varies as 1/co 5 . Thus we might expect, for a given 
experimental arrangement, very large errors at low 
frequencies and negligible errors for sufficiently high 
frequencies. 

These conclusions can be checked against the 
results obtained by Wilson [1951] whose report contains 
all of the necessary details. The circles in figure 1 
show a ob Ja, the ratio of the observed to the known 
attenuation in water 3 at 4 °C, as measured in a 12- 
liter borosilicate boiling flask. The errors at the 
low-order modes are seen to be enormous; the observed 
attenuation is too high by a factor of more than 100 
at 16 kHz (mode 3). Even at the ninth mode, about 46 
kHz, the factor is about 4. At the higher-order modes 
the factor remains constant at about 1.4. 

Wilson's report [1951] gives all the data needed to 
calculate the error from eq (33). There is some question 
about the Q of borosilicate however. The value most 
often cited, 1200, is taken from Mason [1950]. Wilson 
[1951] measured it as 2200 for longitudinal vibrations. 
Measurements made in our laboratory on two spherical 
flasks vibrating in a vacuum in the breathing mode 
give about 2700 at 4 °C. 

The full lines in figure 1 show 1+ (D e /Dl) from eq 
(33) for (Vs of 1000 and 3000. It is clear that eq (33) 
explains the low-frequency errors fairly well at least 



Equating these gives as the condition of resonance 



3 The "known values" are those extrapolated from high-frequency measurements on the 
supposition that a/oj 2 is independent of frequency. 



29 



100 



50 



20 



o 

z 
< 



10 



1 1 1 1 


1 1 


1 1 1 1 




- 


\ 








— 












- 


\ \° 






- 


- 


w 


oo 

o 
o 




- 




1 1 l^fc 


oo o ° 




@ I 



10 20 50 100 
FREQUENCY, KHZ 



200 



FIGURE 1. Losses in water at 4 °C, 12-liter spherical Pyrex flask. 

The circles show the ratios of the measured to the known values of the attenuation accord- 
ing to Wilson [1951]. The curves show 1+ (De/D,.) from eq (33); the upper curve is for 
(?£ = 1000 and the lower for <? fc = 3000. 

in a qualitative way. It must be remembered that there 
are other sources of error. Some of these will be dis- 
cussed later. The constant error at high frequencies, 
which remains unexplained, is discussed in section 5. 

3.2. Cylinder, Radial Axial Modes 

The calculation is carried out much as for the sphere. 
In the case of pressure-release ends, the result is 

De = 



1 /wo 

c L QE(o(alo) 2 ) \o) 



R 



R>+ {k r aY[l- (Wco) 2 ] 2 

(34) 



In deriving eq (34) the point admittance of the envelope 
was calculated in a manner similar to that for a sphere 
[eq (31)]. This is probably good enough for the (/, 0, 1) 
modes, but dubious for the (/, 0, n>\) modes. A 
more ambitious attack was made by Lambert [1953], 
but his results are not applicable here. The point 
admittance of the liquid, at the interface, is 

Ua = —dip/da _ k r Ji(k r a) 
p a i(opL$a io)p L Jo(k r a) ' 

so that the condition of resonance is 



k r aji(k r a) _ R 

Jo(k r a) 1 — (wo/w) 2 

in which again, R = p L a/p E t. 



(35) 



The main point of eq (34) is that the envelope losses 
for a cylinder behave like those for a sphere, that is, 
they can be very large or very small for a given setup, 
depending on whether the mode numbers are very low 
or very high. 

4. Losses in the Viscous Boundary Layer 

The viscous boundary-layer losses are calculated 
on the assumption that the boundary layer thickness, 
d v = (2/x/p/.co) 1/2 , is small relative to the interior 
dimensions of the resonator. The method is very old; 
a good simple account is given by Lambert [1953]. 
The viscous loss per unit area at the interface is given 
by 

-jf=y/a>PLiJLl2ul (36) 

in which u t is the tangential component of the fluid 
velocity at the interface as calculated from the 
potential. 

It follows from eq (36) that if Ut = everywhere at 
the interface, the boundary-layer losses are zero. This 
is the case for the radial vibrations of an ideal sphere. 
In practice, the filling hole prevents the motion from 
being radially symmetrical, but the resultant field is 
too complicated to work with analytically. 

More generally, the boundary-layer losses vanish 
if the pressure is uniform, and in particular, zero, on 
the interface. For then the tangential component of 
the pressure-gradient, which is proportional to |co^|, 
is zero. It would appear then that other things being 
equal, the boundary-layer losses are smaller the thinner 
and lighter the envelope, that is, the nearer the fluid 
boundary is to pressure release. 

For the radially symmetrical modes of a cylinder 
having pressure release ends, the potential is 



ilj = Jo(k r r) cos — z 



(37) 



in which n is an odd integer. The only nonvanishing 
component of Ut is (u z ) z=a = — dilj/da. This is calcu- 
lated from eq (37) and substituted into eq (36). Inte- 
gration over the lateral area gives 



D B =?j^(a>p,ixl2)^JUkra). 



(38) 



As expected, D v = if the pressure, proportional to 
Jo(k r r), is zero at r=a. We compare this to the loss 
in the liquid itself, which according to eq (14) is given 
with sufficient accuracy by 



(A + 2/*) A; 4 i//W 



and get 



{• 



Db 



3n 2 TT 2 (2cop,,/^) 1 / 2 Oci {k r a) 2 [l-(coolw) 2 ] 
4(2L) 2 k*a ' a R 2 +(k r a) 2 [l 



(Wco) 2 ] 2 
(39) 
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In obtaining eq (39) we have again made use of the 
resonance equation, eq (35) and written a c \jot for 
4^/3 (X+ 2/x). 

Comparison of eq (39) with eq (34) shows that the 
ratio of the boundary-layer loss to the anelastic loss 
in the envelope varies as o)^ 2 . Thus we might expect 
that the loss in the envelope predominates at low, and 
that in the boundary layer at high, frequencies, and 
that both are eventually negligible. 

Computation shows that in a practical case the mode 
number may be high enough so that although the loss 
in the envelope is negligible, the viscous boundary 
layer loss is still considerable. This viscous loss might 
then be calculated by the method outlined above and 
the results corrected accordingly. Mulders [1950| made 
some interesting measurements on water in a cylinder 
28 cm O.D. by 10 cm high. The reverberation method 
(excitation over a band of frequencies) was used at 
mean frequencies of 750, 1010, and 1500 kHz, corre- 
sponding to very high-order, densely packed modes. 
Mulders calculated the viscous boundary loss on an 
"equipartion" basis for a rigid wall and obtained good 
results even at the lowest frequency of 750 kHz, for 
which the boundary layer loss was almost half the total. 



5. Discussion 

The resonator-decay method for the measurement 
of sound absorption in liquids is known from experience 
to have serious drawbacks. Our analysis, rough though 
it is, elucidates the mechanism of some of the losses 
hitherto considered adventitious, and thus provides a 
rationale for the design of experiments. 

We have shown that the basic theory of the method, 
usually taken for granted, is in fact only approximate 
for realistic boundary conditions, but that the errors 
resulting from its application are small for sufficiently 
high mode numbers. Much more serious are the ane- 
lastic losses in the envelope. These are relatively enor- 
mous for low-order modes, but become negligible for 
modes of sufficiently high order. The same is true for 
losses in the viscous boundary layer. Although we can- 
not calculate these in the case of a sphere, for which 
the boundary condition is not uniform owing to the 
presence of the filling hole, we get a fair idea of the na- 
ture of the effect from the analysis of a cylinder. In any 
case, this loss can be minimized by minimization of the 
mass of the envelope relative to that of the liquid. 

In our example we have chosen water as a test liquid. 
Water, having relatively low values of a and of a/a c i, 
gives a "worst-case" type of example, and the unwanted 
losses are relatively much less in complex liquids hav- 
ing high viscosities and high ratios of compressional to 
shear viscosity. The envelope material used in the ex- 
ample, Pyrex, having a Q in the range 2000 to 3000, is 
typical. A higher Q is of course better, but for a given 
setup, a large increase in this Q results in a disappoint- 
ingly small reduction in the lowest frequency at which 
the losses are acceptable. 

The major problems of the method are inherent. The 



experiments are easier to perform at the lower-order 
modes because the modes are well separated and easier 
to identify, and because the demands on the accuracy 
of the geometry are tolerable, but the lower the mode 
number, the greater the unwanted losses. Even when 
the method is used differentially, using a "standard" 
liquid for which the losses are supposedly known, most 
workers have had to avoid the lower-order modes. 

Another matter of concern is that in the experiment 
we have used as an example [Wilson, 1951], there is a 
residual loss of about 40 percent at the high-order 
modes. It is doubtful that this is due to conduction 
of sound by the supports and the discrepancy remains 
unexplained. Moen [1951] in an experiment similar 
to Wilson's (water, 12-/ spherical flask) also found a 
residual loss of about 40 percent. However, he found 
greater losses in smaller flasks, and by plotting the 
decay rate against reciprocal radius and extrapolating 
the resulting straight line to l/a = 0, obtained values of 
alp in agreement with those measured in the tens of 
megahertz range even at 150 kHz, the lowest frequency 
he worked with. Similar results, namely that the 
residual loss is proportional to the surface-to-volume 
"ratio" of the resonator, were obtained by Moen for 
several other liquids. Thus it appears that this loss is a 
boundary effect not accounted for by our analysis. 

Mulders' [1950) results, described in section 4, are 
of great interest. He considered that the good accuracy 
of his calculations of viscous boundary layer loss, 
based partly on the supposition of a rigid wall, furnished 
empirical proof that the walls behaved as though rigid, 
and called this an "amazing fact". According to our 
analysis (sec. 3.) this behavior is to be expected under 
the conditions of Mulders' experiment. 4 

The lively interest in possible laboratory measure- 
ments of absorption at low frequencies, say down to 1 
kHz, has recently intensified. The results of our 
analysis make us pessimistic about this possibility, 
especially for low loss liquids, except perhaps by 
methods as yet undreamt of. Mulders' [1950] experi- 
ment, properly scaled to 1 kHz, would require a tank 
having dimensions of tens of meters. To keep such a 
volume clean, free of air bubbles, and at uniform 
temperature would tax the resources of even a labora- 
tory which could afford to build the tank. 



The (^-values of several empty flasks vibrating in the 
breathing mode were measured by Carl E. Tschiegg. 
He also made numerous measurements, not extensive 
enough to report here, on the vibrations of a liquid 
filled barium-titanate cylinder. These served to verify 
the general nature of the variation of the losses with 
mode number as developed in the text. 

The work was supported in part by the Office of 
Naval Research. 



4 However, we have considered only axisymmetric modes. An unsymmetrical mode of 
the liquid, such as occurs in Mulders' experiment, drives a flexural mode of the envelope 
at a frequency much higher than its natural frequency. Hence its impedance is very high 
and its motion very low. Of course this is only an informed conjecture. 
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